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Abstract —This paper thoroughly explores the fundamental 
Interactions between cell association, cell load and throughput 
in a green (energy-efficient) small cell network in which all base 
stations form a homogeneous Poisson point process (PPP) of 
intensity Ab and all users form another independent PPP of 
intensity Ay- Cell voidness, usually disregarded due to rarity in 
cellular network modeling, is first theoretically analyzed under 
generalized (channel-aware) cell association (GCA). We show that 
the void cell probability cannot be neglected any more since it 
is bounded above by exp(—A U( /A b) that is typically not small in 
a small cell network. The accurate expression of the void cell 
probability for GCA is characterized and it is used to derive the 
average cell and user throughputs. We learn that cell association 
and cell load Au/Ab significantly affect these two throughputs. 
According to the average cell and user throughputs, the green 
cell and user throughputs are defined respectively to reflect 
whether the energy of a base station is efficiently used to transmit 
information or not. In order to achieve satisfactory throughput 
with certain level of greenness, cell load should be properly 
determined. We present the theoretical solutions of the optimal 
cell loads that maximize the green cell and user throughputs, 
respectively, and verify their correctness by simulation. 

Index Terms —Green Communication, Small Cell Network, 
Cell Association, Cell Load, Throughput, Stochastic Geometry. 

I. Introduction 

In recent years, we have witnessed a trend that mobile de¬ 
vices, such as smart phones and tablets, have been proliferating 
and relentlessly penetrating our daily life. Such powerful hand¬ 
sets have created a new dimension of transmitting information 
over newly developed wireless technologies. Accordingly, 
this situation has urged customers to expect more and more 
throughput to suffice their versatile wireless demands. To boost 
network throughput, one of the most effective approaches is to 
increase the spatial frequency reuse by deploying more cells 
in a given area so that there are fewer users sharing time 
and/or frequency resources in a cell |T]. Furthermore, as cell 
intensity (density) increases, users are able to connect to a 
closer base station and thus desired signals suffer less path 
loss. Therefore, a small cell network consisting of highly dense 
picocells and/or femtocells has gained much attention since 
small cell base stations (BSs) can be deployed indefinitely 
and the achievable throughput seems not to have a hard limit. 

Migrating the topology of traditional macro cells to the 
one of small cells is the trend for the next generation of 
cellular systems since a small cell BS has much less power 
consumption than a traditional macro cell BS whose power 
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consumption typically accounts for 80% of the total power 
consumption of a cellular system Cellular networks with 
a small cell topology absolutely help the entire information 
and communication technology (ICT) industry lessen CO 2 
emission, which is estimated to account for about 6% of the 
global emission in 2020 0- Small cell BSs also have an 
advantage of low operational cost and hardware complexity. 
For example, they can quickly switch their different power 
control modes for saving power 0. Due to dense deployment, 
however, small cell networks seriously suffer the void cell 
issue that is hardly noticed in the traditional marco BS 
network. The void cell problem originates user-centric cell 
associatioif] and it could give rise to no users in a cell. To 
explain this cell voidness phenomenon, consider there are 
mxn users uniformly distributed in a unit-area network which 
is tessellated by n equal-sized cells of BSs and thus the user 
intensity is mxn and the BS intesnity is n. The probability of 
having no users in a cell is (l — -) which is close to e -m 
as n is very large. This tells us that the void cell probability is 
not negligible any more in a small cell network because such 
a network usually has a small m, i.e., the ratio of the user 
intensity to the BS intensity is not large. 


A. Motivation and Prior Work 

Without modeling the impact of void cells, the analyses 
of network performance metrics, e.g., coverage probability 
and average throughput, are apparently underestimated since 
the void BSs do not contribute any interference. Despite the 
important phenomenon of void cells, there are very few works 
that study and model the effect of void cells. Almost all the 
prior works on cell association in cellular networks overlook 
this problem 0-(TO']. References GD-Q3 do consider the 
void cell impact on their models, however, the void cell 
probability they used is only valid for nearest BS association. 
Although more recent works in 1141 and ] 1 5) also incorporate 
the void cell issue in their models, their cell void probabilities 
are derived based on the cell association scheme with constant 
biased weights and they are no longer valid for the cell 
association schemes with random biased weights. 

Since cell voidness essentially cannot be completely avoided 
and small cell BSs can easily switch their power control 
modes, the void BSs can be put in a “dormant” mode for 


'User-centric cell association means that every user in the network tries to 
associate with its best service BS via some cell association (or called user 
association in some literatures) schemes. For example, if all users are looking 
for a BS that is able to provide the long-term strongest signal power, they 
will associate with their nearest BS if fading and/or shadowing effects are 
completely removed from their channels. 
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saving power. Such energy-saving strategy entails another 
question, that is, how much denseness of BS deployment a 
small cell network needs to optimally exploit the performance 
of the energy-saving strategy. A very much dense small cell 
network is absolutely not green (energy-efficient) because its 
per-unit power throughput is significantly weakened due to too 
much interference. Previous works on designing the BS inten¬ 
sity with greenness in a heterogeneous or small cell network 
are still fairly minimal. The work in | fl6) numerically showed 
that there exists an optimal ratio of picocell to macro BS 
intensity that maximizes energy-efficient throughput without 
giving a rigorous theoretical explanation. In fT7) , a closed- 
form upper bound on the optimal BS intensity subject to 
the constraints on the user outage rate was found without 
considering the energy consumption of a BS. 

The two aforementioned works did not study the optimal BS 
intensity from the viewpoint of an achievable green through¬ 
put. Although the average achievable throughput of a link was 
investigated in |5| and G3> it fails to characterize how energy 
is efficiently used to transported information on the per cell 
or per user basis. References 0 and GD also present some 
approximated average per-user achievable throughputs in a cell 
that do not consider the void cell issue and are only valid for 
nearest BS association with constant biased weights. In J20) , 
an energy-efficient throughput metric was proposed based on 
the concept of outage capacity. None of these prior works 
on the average per-cell or per-user throughputs investigate the 
optimal cell load problem as well as consider the void cell 
issue in their models. 


whereas the average user throughput is a per-user average 
throughput metric in a cell. The near closed-form expressions 
for these two throughputs parameterized by cell load Au/Ab 
are derived and they are the neatest results than any other 
similar throughput results in the literature. Most importantly, 
they directly reflect that the GCA scheme favoring large 
channel power can significantly benefit them. Also, we show 
that there exists a unique optimal cell load that maximizes 
the average user throughput and it is accurately equal to the 
fixed point of a special function pertaining to the void cell 
probability. 

A simple green power control, which switches a BS between 
active and dormant modes depending on whether the BS is 
void or not, is used to help the network save energy. The 
green cell and user throughputs are proposed and they aim to 
characterize how green the per-cell and per-user throughputs 
can be in the network. We show that both the green cell and 
user throughputs can be maximized by a sole optimal cell 
load that can be accurately acquired by calculating the fixed 
point of the derived special function. Once the optimal cell 
load is found, the optimal BS intensity corresponding to any 
given user intensity is also readily obtained and deploying 
BSs based on this intensity can attain the maximum energy 
efficiency of information transportation. From these derived 
results, we can conclude that green power control schemes that 
keep large power consumption difference between active and 
dormant modes and cell association schemes that favor high 
channel power both significantly improve the optimal cell load 
and green throughputs. 


B. Contributions 

To characterize the void cell probability in a small cell 
network with spatial randomness, in this paper we consider 
a single-tier small cell network in which all BSs form a 
homogeneous Poisson point process (PPP) of intensity Ab 
and all users form another independent PPP of intensity 
Ay. The fundamental interaction between user-centric cell 
association and void cell probability is first delved. All users 
associate with their serving BS by the proposed generalized 
cell association (GCA) scheme that is able to cover several 
cell association schemes, such as nearest cell association, 
maximum received power association, etc. We theoretically 
show that the achievable lower bound on the void probability 
of a cell is exp where Au/Ab is termed “cell load”. An 

accurate closed-form expression for the void cell probability 
of GCA is derived under a more practical channel model that 
incorporates the effects of path loss, Rayleigh fading and log¬ 
normal shadowing. The derived void cell probability provides 
us some insights on how the void cell probability is affected 
by channel impairments and how its theoretical lower bound 
can be achieved by GCA. To the best of our knowledge, the 
fundamental relationship between random cell association and 
void cell probability is first studied in this work. 

Under the assumption that all the associated BSs form a 
thinning homogeneous PPP, the average cell and user through¬ 
puts are proposed to capture the impact of cell voidness. The 
average cell throughout is a per-cell average throughput metric, 


C. Paper Organization 

The rest of this paper is organized as follows. In Section [II] 
the system model of the small cell network with generalized 
cell association is introduced and the preliminary results of 
the random conservation property of a general PPP and the 
power consumption model of a small cell base station are 
presented. Section III studies the void probability of a cell, 
cell load, average cell and user throughputs. The results of 
optimal green cell load and green cell and user throughputs 
are elaborated in Section [TV| Finally, Section [V] concludes our 
important findings and observations in this paper. 


II. System Model and Preliminaries 
A. Network and Cell Association Models 

We consider an infinitely large and planar small cell network 
in which all users form a homogeneous Poisson point process 
(PPP) <f>u °f intensity Ay whereas one-tier (small-sized) base 
stations that provide service to all users independently form 
another marked homogeneous PPP <1>b of intensity Ab given 
by 

$b : Bi G G K ++ ,Ci C R 2 , 

V, G {0,1},V*GN+}, (1) 

where £/ denotes base station i and its location, II, is used 
to characterize the downlink fading and shadowing channel 
power gain from Bi to its service user, the cell region of II, 
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is represented by Ci, Vi is a void cell index that indicates 
whether or not Cj has a user, i.e., whether Ci D $u = 0 is true 
or not - Vi is equal to one if Ci n $u 0, otherwise zero. In 
order to capture stochastic behavior of the channel power in 
the downlink, we assume that all channel power gains Hi s are 
i.i.d. random variables and their probability intensity function 
(pdf) that characterizes the composite effect of Rayleigh fading 
and log-normal shadowing is given by ED 


= 


1 f°° 1 

— , / —=■ exp 

a/ 2 t uj'I J o+ x 2 



(\nx-p s ) 2 \ 
2a 2 ) 


( 2 ) 

where p s and a 2 are the mean and variance of log-normal 
shadowing, respectively. 

Without loss of generality, our following analysis will be 
based on a typical user Uq located at the origin. Each user 
associates with a base station in $b by using a generalized 
(channel-aware) cell association (GCA) scheme, i.e., user Uq 
associates with its serving base station Bq via the following 
scheme 


Bq = arg sup (WjfTj||.Bj|| a ) , (3) 

Bi G<t>B 


where all Wf s are the i.i.d. random association weights for 
BSs, a > 2 is the path loss exponent and ||.Bj|| denotes the 
Euclidean distance between Bi and the origin. The motivation 
of proposing the GCA scheme in |3]l is two-fold. First, GCA 
can generally cover cell association schemes with deterministic 
and/or random association weights. For example, if the channel 
power gain Hi is available and W, = 1/Hi, GCA reduces 
to nearest BS association. GCA becomes maximum received 
power association provided that all Wf s are the same constant. 
Second, GCA can be viewed as an adaptable cell association 
scheme and it is suitable no matter whether BSs can acquire 
the mean channel power gains of users in time or not. For 
example, BSs may not be able to estimate the mean received 
power from non-stationary users that are moving very fast ED- 
GCA is essentially a “user-centric” scheme, that is, it is able to 
ensure every user to connect to a certain base station. No users 
are blocked out of the network. Nonetheless, user-centric cell 
association cannot, as we will show later, guarantee every cell 
is associated with at least one user, i.e., the probability of a 
void cell (a cell without users) is always bounded above zero 
and non-negligible, especially in a small cell network with 
a large BS intensity. A simulation example of the void cell 
phenomenon in a small cell network with Voronoi tessellation 
is illustrated in Fig. [D 


B. Random Conseri’ation Property of a General PPP 

In this subsection, we introduce the random conservation 
property of a general PPP, which specifies how the intensity 
measure of a PPP is changed after all points of the PPP are 
transformed by i.i.d. random mapping matrices. This property 
is elaborated in the following theorem. 


Theorem 1 (Random Conservation Property of a General 
PPP). Suppose •V' is a marked general PPP of intensity 
measure A' on R d , which can be expressed as follows 


= {{X', Tj) : X[ G R d ', T t G R 


,V*GN+}, (4) 



Fig. 1. An illustration example of void cells in a small cell network. BSs 
(blue crosses) and users (brown dots) are two independent PPPs and the cells 
are created by Voronoi tessellation. The intensity ratio of users to base stations 
is Au = 2. 


where X[ denotes node i and its location and Tj : K d —> K d 
is the non-singular mapping matrix (operator) of node X[. For 
all i j, Tj and Tj are two different random matrices and 
their corresponding elements are i.i.d. random variables. Let 
<3E> be the mapped point process on R d generated by using the 
random mapping matrix of each node in <!>', i.e., it is defined 
as 


4 = {Xi 4 Tj(X') : X[ G <f> / , Tj G K dxd ',Vz G N+}. (5) 

For any d-dimensional bounded Borel set A C R d and d 7 - 
dimensional bounded Borel set A' C R d , ifv d '(A') = Vd(A), 
then <t> is a general PPP of intensity measure A given by 

1 


A(.A) = A'(.4')E 


( 6 ) 


Vdet (TTT) 

where T 1 is the transpose of T. If 4> 7 is homogeneous with 
intensity A 7 , •!> is also homogeneous and has the following 
intensity 

1 


A = A'E 


(7) 


ydet(TTT) 

Proof: See Appendix [A] ■ 

Remark 1. Theorem 7] is a generalization of the consen’ation 
property in [22^ [23], In a special case of d = d 7 = 2, all 
points in $ are mapped from their corresponding points in 
a homogeneous PPP 4> 7 by scaling them with i.i.d. diagonal 
random matrices Tj = diag(Ti,Ti). In this case, A is equal 
to A'E[T -2 ]. 


The random conservation property can significantly reduce 
the complexity of analyzing the statistics of some performance 
metrics induced by a PPP, especially a homogeneous PPP with 
i.i.d. marks. Hereupon the GCA scheme in ([3]» can be rewritten 
as Bq = arginf^G$ B \\{WiHi)~^Bf\\ which can be further 
simplihed by Theorem [D as 

\\B*\\±(WH)^ (.inf ||All), 


dxd 7 


( 8 ) 
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where = stands for equivalence in distribution, Bi = 
( WjH ,)“ o Bj and <E>b is a homogeneous PPP of intensity 


ArE 


{WHY 


Hence, the distri- 


based on ([7} in Theorem 
bution of the distance between user Uq ana its serving BS Bq 
can be instead equivalently found by the distribution of the 
{WH)& -weighted distance from the origin to the nearest BS 
in the new PPP <f>B. In other words, the random transformation 
property can transform GCA into another form of nearest BS 
association , which significantly simplifies the analysis of the 
performance metric of GCA, such as coverage/outage prob¬ 
ability, since many existing results of nearest BS association 
can be applied in this context by simply modifying them with 
an updated intensity of the BSs. 


C. Power Consumption Model of a Small Cell BS 

Although the cell voidness issue mentioned in Section II-A 
seems inevitable, its impact on network energy consumption 
can be alleviated by applying green power control at BSs. 
The green power control scheme for a small cell BS has two 
modes - the active mode is for non-void BSs, whereas the 
dormant mode is for void BSs. For base station Bi, the power 
consumption of such green power control can be characterized 
by the following expression | |24) 

(9) 


= Won + SP t ) + (1 - Vi)P oS , 


where (H, is the power consumed by BS /i,, P on is the power 
consumed by the hardware of a BS in the active mode, P t 
denotes the constant transmit power of a base station, 6 > 0 
is a scaling constant for the transmit power which usually 
depends on the power amplifier used by a BS, and P 0 g denotes 
the power consumption of a dormant BS. Also, we assume 
P on > P 0 g since it is usually the case in practice. The average 
power consumption in |9| can be respectively written as 

E[*]=(l -p®){Pon+6Pt)+p 9 Poff GO) 


in which p$ = P[Vj = 0] represents the void probability of 
a cell and its analytical result will be elaborated in Section 


III 


Note that is certainly affected by cell association 
schemes adopted by users. In Section [ill} we will investigate 
how to characterize p@ for GCA in ([7|i from a fundamental 
connectivity perspective. In addition, the transmit power P t 
of all BSs is assumed to be able to provide the minimum 
received power P m j n to their associated users by compensating 
the mean path loss, i.e., P t = P m i n inf{E[||Bo 11“]}- According 


to 


| Pq|| in <[8]>, we have 


nmn =e 


=E 


2n\ B E[{WH)«]{WHy 


r l+a e -»r X B E[(WH) 5 ] ( WH) - 


dr 


r(i + f) 


{■K\ B {WH)-il&[{WH)i})f 


> 


r(i + f) 

(7rA B C) f ’ 


where T(x) = L°° t x 1 e 4 d t is the Gamma function and £ = 


E 


(WH)i 


E 


(WH)- 


Pt = 


> 1. Thus, P t is given by 


PminlXl + f ) 


(ttAbC)- 


( 11 ) 


Note that P t in (11) is implicitly assumed not to exceed the 
maximum transmit power provided by a small cell BS for a 
given P min and A B . By plugging © into the explicit 
result of E[\f r ] is acquired, which will be used to define the 


green cell and user throughputs in Section IV 


III. Void Cell Probability, Cell Load and Average 
Throughputs 


A. The Void Probability of a Cell for GCA 

As pointed out in the previous section, the GCA scheme 
or any “user-centric” association scheme cannot guarantee 
that there are at least one user in each cell, i.e., a void 
cell could exist in the network. This phenomenon can be 
intuitively interpreted by using a Poisson-Dirichlet (Voronoi) 
tessellation for a PPP of BSs. Suppose the Voronoi tessellation 
is used to determine the cell of each BS in $ B and all users 
adopt the nearest BS association scheme to connect with their 
serving BS. This nearest associating process can be viewed 
as the process of dropping all users in $u on the large 
plane consisting of the Voronoi-tessellated cells formed by $b- 
Under this circumstance, the probability mass function (pmf) 
of the number of users in a cell of $ B can be expressed as 


p n ^P[$u(C)=n]=E 


(Au u{C)r r _ x ^o 

I 


n\ 


( 12 ) 


where C denotes a Voronoi cell of a BS in <I> B , ‘by (C) 
represents the number of users in cell C, and v{C) is the 
Lebesgue measure of C. 

Unfortunately, the theoretical result of p n in © is un¬ 
known since the pdf of a Voronoi cell area is still an open 
problem |22) . However, it can be accurately approximated by 
using a Gamma distribution with some particular parameters 
(|22|, |25| . Reference | [25) suggests the following Gamma 
distribution for / y (c) (x)- 


U{C){x) 


{p\zx)P , XrX 
T{p)x 


(13) 


and p = | can achieve an accurate pdf of a Voronoi cell area. 
Substituting (fT3| into (fl2|) yields the following result: 


Pn 


A U (Mb) P f°° n+p-1 e -(pA B +Au)x d 

n! T{p) Jo 

1 T(n + p) / Au/Ab \ n f Au 
n\ V{p) \p + \ u/A b / V Mb/ 


(14) 


We call the term Au/Ab in ( fT4] > the cell load of the network 
since it represents the mean number of users in a Voronori- 
tessellated cel Q] Hence, the pmf of the number of users in 
a cell for nearest BS association has an accurate closed-form 
expression. The void probability of a cell, can be found 
by Pn with the case of n = 0, which is the last term of p n in 
( fl4| ), i.e., 

m 4p [K = 0]=(i + Ll) " (15, 


2 Although the pdf of Voronori cells for a homoge neo us PPP of inten¬ 
sity A is unknown, its mean can be shown as 1/A |22| . Thus, cell load 
= E[AuKC0] = Au/Ab. 
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and this indicates that the intensity of the void BSs is Agpg 
that is not negligible especially when cell load Au/Ab is 
small. Most importantly, the results in © and © are no 
longer accurate for all non-nearest BS association schemes 
since users do not necessarily associate with their nearest BS, 
such as the GCA scheme proposed in Section [77] However, an 
accurate void probability of a BS for GCA can be derived 
as shown in the following. By Jensen’s inequality, the lower 
bound on pg in (fl5| is given by 


PH = E 


e -Au^(C) 


> exp 



(16) 


This lower bound on the void probability reveals three crucial 
implications: (i) the void probability of a cell is always 
bounded above zero such that there could be a certain num¬ 
ber of void BSs in the small cell network; (ii) nearest BS 
association cannot achieve this lower bound since its void 
probability in © is always larger than the lower bound; (iii) 
from an energy-saving perspective, the lower bound represents 
the minimum percentage of void BSs in the network that can 
be turned off to save energy. Later, we will theoretically show 
that, this lower hound can be achieved by GCA. To derive an 
accurate void probability of a cell for the GCA scheme, we 
approach this problem from a fundamental connectivity point 
of view and derive the bounds on pg as shown in the following 
theorem. 


Theorem 2. If all users adopt the GCA scheme defined in ([3]) 
to associate with a BS in $b defined in 0. the bounds on the 
void probability of a cell are given by 

( 1 + ^?) <17) 



Fig. 2. The void probabilities of a cell for nearest BS association (Wi = 
1 /Hi, p = and maximum received power association (Wi = 1, p > /). 
The network parameters for simulation are path loss exponent a = 3.76, 
shadowing mean /r s = 0 dB, Ay = 370 users/km 2 . 


This indicates that the void probability of a BS is reduced 
when more users join the network under a given BS intensity. 
In other words, when the network has a large user population 
the efficacy of reducing pg by using large ( is apparently 
undermined such that the performance of GCA is similar to 
that of nearest BS association in this case. 

Although the bounds on pg are characterized, an accurate 
result of pg is still needed since it will help us understand how 
many BSs per unit area are void and they should be switched 
to the dormant mode. The following proposition renders an 
accurate heuristic result of pg for GCA. 


and ( = E 
Section I//- Cl 


{WHY 


E 


(WHf 


> 1 as defined in 


Proof: See Appendix |B| ■ 

According to the proof of Theorem [2] the lower bound on 
pg is derived while considering the completely independence 
exists between the non-associated events of a BS whereas the 
upper bound is obtained by approaching the opposite case, i.e., 
all non-associated events of a BS are caused by considering 
all users in a cell located the same farthest distance from 
the BS and thus they are highly correlated. Accordingly, it is 
reasonably to conjecture that the upper bound is tightly close to 
the lower bound provided that the cross-correlations between 
all non-associated events are significantly weakened. On the 

other hand, mathematically we know lim^oo ^1 + = 

exp anc * t ^ lus t ^ le bounds in ( fl7] > are fairly tight as 

f becomes large. This intuitively reveals that large £ will 
“decorrelate” all non-associated events, which is an important 
observation since it lets us realize that the lower bound on 
pg is easily achieved by enlarging the --fractional moment of 
WiHi. For example, the lower bound on pg can be achieved 
by the maximum received power association scheme provided 
that channels have large shadowing power. In addition, the 
bounds in (fI7| can be tight for the case of a large cell load. 


Proposition 1. The void probability of a BS can be accurately 
approximated by 


p »=( i+ ^r <i8) 

where p = ^E[(W H)i]E[(W H)~ i] = f£ if the GCA 
scheme is used in the network. 

Proof: Since (1 + x)~ l > (1 + x/a)~ a for a > 1, letting 
x = Au/Ab and a = p/C leads to 


1 + 


Au A C 

abc; 


> 



Ab p) 


since p > ( > 1. Based on the accurate approximation of the 
pdf of a Voronoi cell area suggested in [251, pg with p = | 
in ( p~5j ) is an accurate void probability of a BS for nearest BS 
association. Therefore, we can conclude thatpg with p = is 
accurately the void probability of a BS for the GCA scheme 
since such p reduces to f as GCA reduces to nearest BS 
association (i.e., W. L = 1/Hi for all i). ■ 

Fig. [2] shows the simulation result of the void probabilities 
for nearest BS association and maximum received power 
association. For nearest BS association, we have p = | since 
Wi = 1 /Hi, while p for maximum received power association 
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is given by 

- . ^ , exp(^) , (19) 

asm(27r/aj \ cr J 

which only depends on the shadowing power (variance) and 
path loss exponent. First, we see that the void probability 
of nearest BS association is no longer accurate if maximum 
received power association is used. For example, when ^ « 
2, the void probabilities for the nearest BS and maximum 
received power association schemes with 8-dB shadowing 
are 0.2 and 0.14, respectively. The lower bound on the void 
probability is around 0.135. The void cell probability given in 
m indeed accurately coincides with the simulated result and 
is much closer to the lower bound. Next, as expected, large 
shadowing power (i.e., p » 1) indeed makes the void cell 
probability approach to its lower bound. Thus, GCA is able 
to achieve this lower bound. 




H< 


E 


H~ 


B. Average Cell Throughput 

Since void cells could exist in the network and they do 
not contribute network interference, the throughputs obtained 
in the previous works overlooking the void cell problem are 
certainly underestimated if interference is treated as noise. 
To effectively quantify how the void probability impacts the 
downlink throughput of a BS, we define the following average 
cell throughput assuming the channel input is Gaussian and 
interference is treated as noise: 


7c = sup 
b5g$ b 






where Iq = Yh B g<s> b \b* ViHi\\Bi\\~ a is the interference at 
user Uq assuming all spectrum is fully reused and z/(Gq) is the 
Lebesgue measure of cell Cq associated by user Uq. Average 
cell throughput 7c characterizes the trade-off relationship 
between average downlink channel rate and cell size due to 
the void probability - large void probability benefits average 
channel rate owing to interference reduction, yet it increases 
the average cell size of the associated BSs that cover the whole 
network such that more path loss is induced. The explicit result 
of 7c is given in the following proposition. 


Proposition 2. If every user associates with its cell by GCA 
and all the associated BSs are modeled as a homogeneous 
PPP and let Cz(s) = K[e~ sZ ] be the Laplace functional of 
random variable Z, then 


£||b;||“J 0 (s) = E 


G 2 


_G 2 + (l-p 0 )f(s,G 2 ) 
where p 0 is given in ( p~8] >, G 2 = (WH)~iE[(WH) 


£(s, <f>) =S“ < E 


Hi 


r 11 — I - 

a. 



[l -C H (t 


( 21 ) 


and £||b 0 *ii/ 0 ( s ) < i+ti-p^o.o/c where £ = E t G2 ]- For a 

sufficiently large positive integer n, the average cell through¬ 
put in ( |20[ l can be accurately found by 

= A B (l-p 0 ) 

/C (In 2)(1 — 1/p) 

n /*0O i 

/ uJids 

f^Jo [s + e-(v^«r.x 1+ M.)] [1 + (1 - pv)£(s, C)/C] ’ 

( 22 ) 


where u>i = — Ax )] 2 ’ Xi s are ^ le roots °f dle physicists’ 

version of the Hermite polynomial H n (a:) [26], If there is no 
shadowing in the channels, ( |22| ) reduces to 

A B (1—p 0 ) f°° ds 


7c = 


(ln2)(l - 1/p) 7 0 (s + 1)[1 + (1 -p 0 )^(s,C)/C] 

(23) 


and f(s, </) = si 


2tt 


a sin( 27 r/o:) Jo 


-ft 




dt 

1+t^ 


Proof: See Appendix |C| 


Remark 2. In order to make the derivations of 7c tractable, 
we assume the associated BSs still form a homogeneous PPP 
in this paper even though in theory they are no longer a 
PPP due to the correlation between them (see the proof of 
Theorem [2]). Nonetheless, our previous work in showed 
that the associated BS can be still accurately approximated 
by a homogeneous PPP as long as the cell load is not pretty 
small. 


Remark 3. If the coverage probability is defined as the 
probability that the SIR of the user is no less than a threshold 
s £ R++, it follows that £|| B *||cj 0 (s) in ( |2 1 1 ) is the user 
coverage probability in a cell. The result in © is much neat 
and more general if compared with the coverage probability in 
prior works since it indicates how void cell probability, fading 
and shadowing affect user’s coverage. 


The average cell throughput in ( |22| > is the nearest closed- 
form result than any other similar average throughput results in 
literature (such as @ ©)■ Its physical meaning is the per-cell 
spectrum efficiency of quantifying how the average spectrum 
efficiency of associated BSs changes along their average 
cellular size 1/A B (1 p 0 ) under the GCA scheme. Especially, 
it is the first cell-wise throughput result including the impact 
of the void cell probability. It possesses several important 
meanings that are addressed in the following, respectively. 

. BS intensity dominates 7c in heavy cell load: For the 
case of heavy cell load (i.e., Au/A B 1), we have p 0 « 
0. Intuitively, decreasing p 0 increases more interference, 
however, it also induces less path loss in desired signals 
as well as the average cell size of the associated BSs. In 
this context, 7c is approximated by 


Tc 


Ab 


(ln2)(l — 1/p) 

n pOO 

E 


w,:ds 


o (s + e-(v / 2^^+^))(l + ^(s,C)/C) : 


(24) 


which indicates 7c becomes a linearly increasing function 
of A b such that the asymptotic scaling law 7c G 0(A B ) 


























7 


holds. Thus, under the heavy cell load situation deploying 
more BSs can directly increase 7c- 
. User intensity dominates 7c in light cell load: For the 

case of light cell load (i.e., Au/Ab -C 1), we have ss 1 
and 7c is approximated by 

rc “(ln2)(l-l /p) X 

A r°° _ Agt^ds _ 

hi ( s + -))(Ab + Au £(s, 0/0 ’ 

(25) 

For a given Ay, we can observe that 7c is a monotonic 
increasing and concave function of Ab and thus 7c £ 
0(Ay) as Ab — > oo. For a given Ab and Ay —f 0, we 
also have 7c £ 0(Ay). Accordingly, making more users 
join the network can efficiently boost 7c in the light cell 
load situation. This scaling law corrects the flawed scaling 
result of the average throughput in the previous work 
without considering void probabi I i t\j^] 

• Cell association significantly impacts 7c: Cell associa¬ 
tion schemes have a considerable impact on 7c since pg 
and £(s, Q in ( |2T| both are affected by cell association 
weights. Specifically, a cell association scheme with fairly 
large £ certainly benefits 7c since decreases 

considerably and (1 — pg) almost remains a constant 
in this case. General speaking, nearest cell association 
(WiHi = 1 and ( = 1) is unable to achieve higher 7c 
than any other cell association schemes (with random 
WiHi and ( 2> 1). In maximum received power asso¬ 
ciation, for example, cell association results are highly 
impacted by ( that exponentially increases with shad¬ 
owing power as shown in © such that 7c becomes 
much larger in a higher shadowing power environment. 
In this context, cell association schemes that favor BSs 
with stronger channel power surely benefits 7c- 
A numerical illustration of 7c for nearest BS association 
and maximum received power association is shown in Fig. [3] 
As can be seen, the analytical results in ( |23| with n = 6 per¬ 
fectly coincide with their corresponding simulation results in 
different cell association schemes and shadowing conditions. 
All simulation curves show that 7c is a monotonic decreasing 
and convex function of the cell load for a given Ay and 
eventually it will approaches to a constant regarding that A|f] 
Most importantly, they indicate that maximum received power 
association evidently outperforms nearest BS association in 
terms of the average cell throughput since it has a larger £, as 
we have already pointed out in above. 

C. Average User Throughput 

Although the average cell throughput found in Proposition 
[2] presents the per-cell spectrum efficiency of the small cell 

3 Using the average rate without considering void probability in [ 51 to 
calculate 7c defined here, only the scaling law 7c S <—> (Ab ) for all 
Ay S ®++ can be concluded. This is certainly not true for light cell load 
since a large portion of BSs are not connected in this case so that they do 
not contribute any throughput in the network. 

4 Note that 7c will increase along the cell load and finally converge to a 
constant if the BS intensity is fixed. This is the desired situation we wish to 
have since 7c and cell load both increase at the same time. 



^—Analytical Tc in (22) (NBA, No Shadowing) 
—©—Simulated Tc (NBA, No Shadowing) 

-Analytical Tc in (23) (NBA, Shadowing, a\ = 3dB) 

—©—Simulated Tc (NBA, Shadowing, <7^ = 3dB) 

-Analytical Tc in (23) (NBA, Shadowing, a] = 6dB) 

—$—Simulated Tc (NBA, Shadowing, = 6dB) 

- Analytical Tc in (22) (MRPA, No Shadowing) 

--^"Simulated Tc (MRPA, No Shadowing) 

-Analytical Tc in (23) (MRPA, Shadowing, a] = 3dB) 

•-^-- Simulated Tc (MRPA, Shadowing, a\ = 3dB) 

-Analytical Tc in (23) (MRPA, Shadowing, irj = 6dB) 

■^"Simulated Tc (MRPA, Shadowing, <rj = 6dB) 


\u/\ B (Cell Load) 


Fig. 3. Simulation results of average cell throughput for nearest BS 
association (NBA) and maximum received power association (MRPA) under 
Rayleigh fading and log-normal shadowing. The simulation parameters are 
a. = 3.76, Ay = 370 users/km 2 , fi s = 0 dB and n = 6 for eq. (22) . 


network with cell voidness, it fails to delineate how much 
of the downlink throughput resource in a cell can be fairly 
allocated to users either in the frequency or in the time domain. 
To characterize the average throughput fairly shared by users, 
we define the following average user throughput 


% = sup 


E 


log 2 1 + 


Hn 


\B* 0 \\~ a 


E[$y(C 0 *)] 


(26) 


where $u(Cq) is the number of serviced users in a non-void 
cell. Note that the distribution of <i»y (C g ) is not equal to that of 
$u(Al) for any bounded Borel set A C M 2 with v{A) = J'(Cg) 
due to the cell voidness issue. The following proposition gives 
the explicit result of 7u- 


Proposition 3. If the shadowing effect is considered in chan¬ 
nels, the average user throughput in \26) can be explicitly 
found as 


% = 


Ab(I-Pb) 

Au (hi 2) 


-Ulids 


i=1 ' 


s -(V2a s x i+/ , s )] [1 + (1 - n )t (S, C)/C] ' 

(27) 


If the shadowing effect is not considered, the average user 
throughput is 


Tv = 


ds 


a b (i -p 0 ) 2 r _ 

Au(ln2) 7 0 (s + l)[l + (l-p 0 K(s,<)/<] 


, (28) 


where u>i, Xi, and f(s, f) are the same as given in Proposition 

m 

Proof: According to the pdf of the area of associated BS 
C* in ( |40| , the pmf pf $u(Cq) is 


P[$u(C 0 *)=n]=E 

_ AS(pA B (l - n )y f c 
n!(l-p 0 )r(p) J 0 


(A u ^(Co))” c _ Auy(c . ) 

n > 1 

n\ 



( ,n-l+P e -(A u +p(l-p 0 )AB)x ( j 2 . 



























where 


PH T(n + p) / Au/A b (1 — Pm) \ n 
(1 -pa) r (p)n\ V a u/Ab(1 -pn) + p) 

where p 0 is given in ( | 1 8| . As a result, the mean number of 
users in cell is 


nMc*o)} 


n 

(i -P0)r(p) 

v^r(n + p) / Au/Ab( 1 — p@) A 
( n — !)! \Au/Ab(1 -pm) + p) 
Au 

Ab(1 - P @) 2 


The average user throughput can be calculated as follows 

T _ A b (1 -p$) 2 Tc = (p- 1)(1 -Pg)Tc 9 „ 
/u Au(E[l/i/(Cg)]) P Au 1 } 

Substituting |22| into ( |29[ i results in \21\ and ( |28[ i. ■ 

The average user throughput in ( [27] ) and ( [28] ) is a linear 
scaling of the average cell throughput normalized by the 
user intensity such that it is virtually the function of cell 
load Au/Ab- For the asymptotic behaviors of 7Tj- we have 

7u G 0 (aTO+^OTC)) and 7u 0 in the regime of high 
Au/Ab (i.e., the regime of heavy cell load) since there are 
too many users sharing the throughput resource in a cell. On 
the other hand, in the regime of low Au/Ab (i.e., the regime 
of light cell load), 7u scales like ©(Au/Ab) and 7u —> 0 
since the average area of non-void cells is so large that the 
average throughput shared by users is fairly small. Specifically, 
in moderate cell load there exists a sole optimal cell load that 
maximizes the average user throughput. Knowing this optimal 
cell load is important since it indicates when cell load and 
average user throughput increase at the same time and in this 
case more user throughput can be exploited by deploying fewer 
BSs for a given user intensity. The following proposition gives 
the accurate solution of the optimal cell load for 7u. 


Proposition 4 . The optimal cell load that maximizes the 
average user throughputs in © and ( |28| ) is the fixed point 
of the following real-valued function parameterized by f3 > 1 

P 

, v > 0. (30) 


l(v) = 


1- 1 


-p 


Namely, if the optimal cell load is v*, we have L(v*) = v*. 


Proof: Since the cell load is only contained in the term 
(1 -pa) 2 A B /Au[l + (1 - P0K(s,C)/C]i we its optimal 
value by merely calculating the derivatives of this term with 
respective to the cell load. Let v denote the cell load, V(i>) 
represent this term parameterized by wand pa = 1 — Then 
we have 

V( \ = P0 

W v[l+P 0 *(s, 0/C]’ 

Note that t(s, C)/C is not function of v. Hence, the first-order 
derivative of V(w) w.r.t. v is 

dV _ P0Pa^[2+P0l(s,C)/C] ~Pl\ l +P0l( s >O/C] 
dv v 2 [l +p 0 f(s,C)/C] 2 


_ d P0 _ / -i _l “ 


-p-i 


= p\ +p >o. 


Letting = 0 yields the following differential equation of 


Pa and its solution 


dp0 

dv 


p@[i +P 0 ^(s,C)/C] 


p\ 


= c(p 0 )w, 


w[2+p 0 f(s,C)/C] l+P 0 ^(s,C)/C 

where c is an unknown function of pa that can be determined 
as follows. First, we can assign c(w,pa) = Pa ^/(^ + 
P0^(s,C)/C) for all /3 £ K++ since pa should only depend 

on v, not on f(s, C)/C- Thereby, we acquire pj^ = v. Since 
1 +- 

p'a = -p'a = p 0 p , we know p' 0 G [0, 1] and p 0 < 0. To make 
pj"‘{v) have a fixed point at v*, the first-order and second- 
order derivative constraints on (v), 1 (O)p 0 (O) > 0 

and ~ 1 )P0 //3_2 (P0) 2 + ^0 //3_1 P0 < 0 for a11 v > °, 

must hold. Since they hold only when /3 > 1, L(w) in ( [30] ) is 
obtained. Finally, we can show 


r dV n 
Inn —— > 0 

ov 

~£flPi 


and 


dV n 

lim —— < 0, 

v—too OV 


and the solution of p 0 (w*) — v* = 0 for v* > 0 is unique 
since v* > 0 and p 0 is a positive and one-to-one mapping 
function of v. This follows that L(w*) = (1— p®) 1 ^ = p 0 ' 9 = 
v* and the fixed point of L(w*) is the optimal cell load that 
maximizes 7u. ■ 


Remark 4. For any GCA schemes with p 1, <[30]) reduces to 
L(w) « (1 — exp(— v)) 1 ^ whose fixed point is v = 0. Flence, 
the optimal cell load asymptotically reduces to zero as p goes 
to infinity. 


Cell association also has a pivotal influence on the av¬ 
erage user throughput same as the case in the average cell 
throughput, and this situation can be visually perceived from 
the numerical results of nearest BS and maximum received 
power associations in Fig. [4] Indeed, large shadowing power 
apparently increases 7u as well. There exists only one optimal 
cell load that maximizes 7u as shown in Proposition [4] For 
example, in the case of maximum received power association 
with a = 3.76 and shadowing variance = 4 dB, Proposition 
[4] delivers an optimal cell load « 1 for p = 10.3 and /3 = 2, 
which accurately coincides with the optimal cell load observed 
in the corresponding curve case of Fig. [4] In addition, we also 
observe that the optimal cell load happens later as shadowing 
power gets larger, which further implies that shadowing is 
instructive to the efficiency of deploying cells since higher 
average user throughput is attained by less deployment cost. 


IV. Optimal Green Cell Load and Throughput 

The average cell and user throughputs found in the previous 
section have a common drawback, that is, the energy cost is 
not integrated into them while considering to achieve higher 
throughput by deploying more cells. Deploying more cells 
may not necessarily leads to more total power consumption 
since the transmit power in {H} (that indicates P t oc A B 2 ) 
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Fig. 4. Simulation results of average user throughput for nearest BS 
association (NBA) and maximum received power association (MRPA) under 
Rayleigh fading and log-normal shadowing. The simulation parameters are 
a = 3.76, /i s = 0 dB, Ay = 370 users/km 2 and n = 6 for eq. (27). 


decreases due to less path loss even though other hard¬ 
ware power consumptions increase in this case. In order 
to delve the fundamental relationship between throughput 
and power consumption, we define the following green cell 
throughput Gc (bits/Hz/joule) and green user throughput Q u 
(bits/Hz/joule/user): 


Gc 


Tc 

A b E[*] 


and Q u 


% 

E[\&] 


(31) 


in which Ef'T] is the average power consumption of a BS and 
it is already given in ( fTO] ). The physical meanings of green 
cell and user throughputs are how much throughput can be 
transported by costing unit power consumption. Note that the 
definitions of Gc and Gu are slightly different since Gc aims 
to quantify an energy efficiency on the throughput of a BS 
whereas Gu desires to characterize an energy efficiency on the 
throughput of a user. In this section, we would like to investi¬ 
gate the fundamental relationship between cell load and these 
two green throughputs in order to gain some insights on how 
to deploy BSs with an appropriate intensity for achieving high 
throughput efficiency. The green cell and user throughputs are 
studied in Sections |IV-A| and |IV-B| respectively. 


A. Green Cell Throughput and Its Optimality 

According to ( |3T) and ( |22| ), the green cell throughput can 
be specifically expressed as 


ds 


Gc = 


poo 

Z^i=i Jo [ g+e -(^ s x i +„ o ) ][1+(1 _ p0)<(SjC)/c] 

(In2)(1 - 1/p)[(P 0 „ + SP t ) + p 0 P of f/(l - n )}’ 


(32) 

already given in 0 In heavy cell 
0 and Gc approaches to a constant approximated 


where P t = ~ A" as 
OAbO 2 
load, p 0 ss 0 and y { 

by 


Gc 


E"=i fc 


ds 


o (s+e-(' / 2 <T « a, i+< 1 »))(l-(-^(s,C)/C) 
(In 2)(1 — l/p)(P on + SP t ) 


which can be viewed as the “least” green cell throughput 
merely affected by GCA and BS intensity since the increase 
in the power consumption of the active mode has a direct 
impact on Gc- In this case, deploying more BSs and using 
appropriate cell association schemes with a larger ( can make 
cell throughput greener. In light cell load, p 0 ss 1 and Gc has 
the following approximation 

E n r°° _ds_ 

i=l JO [ s + e -(y/2cr s x i +fj. s )] 

^ ~ (ln2)(l — l/p)P 0 ff ’ 

i.e., Gc = ©(Au/AePoff) as Au/Ab —> 0, and it is the greenest 
cell throughput since any additional power consumption in 
the active mode hardly reduces the cell throughput. However, 
Gc in this case would be very low since it is dominated by 
the light cell load. Lowering the power consumption of the 
dormant mode and increasing cell load can linearly augment 
the green cell throughput. 

The two aforementioned extreme cases indicate that the cell 
load should not be too heavy or too light in order to maintain 
a certain level of cell throughput with satisfactory greenness. 
The exists a unique optimal cell load that maximizes the green 
cell throughput for a given user intensity as shown in the 
following proposition. 

Proposition 5. For a given user intensity Ay, there exists a 
j3 > 1 such that the unique optimal cell load that maximizes 
Gc in ( |32| ) can be accurately found by the fixed point of the 
following real-valued function 


L c (tO =7 tCAu 


( (Pon ~ Ppff) 
^P min T( 1 + 1) 



J2 


(33) 


Namely, the optimal cell load is the unique solution of 
]Lc(t'*) — v* = 0 such that the corresponding optimal BS 
intensity is Ag = Xu/v*. 

Proof: See Appendix [D] ■ 

There are some interesting implications that can be per¬ 
ceived from the fixed point result in ( |33] >. First of all, a nonzero 
optimal cell load does not exist if BSs do not perform green 
power control, i.e., P on = P Q ff, whereas the more power 
is saved (i.e., P on — P Q g is larger), the larger optimal cell 
load and green cell throughput are achieved. Cell association 
schemes with large ( increase the optimal cell load as well 
as the optimal green throughput at the same time. This is 
because large ( reduces the void probability, yet increases 
the coverage probability. Furthermore, cell association has a 
more direct impact on the optimal cell load than green power 
control since the effect of green power control is weakened 
by the power of -. Although green power control and cell 
association both benefit the optimal cell load and green cell 
throughput, green power seems more favorable in practice 
since it is easily implemented and not dependable on the 
communication environments. 
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B. Green User Throughput and Its Optimality 

Substituting ( |27j ) into ( |3Tj ) results in the explicit expression 
of the green user throughput given by 


Q u = 


A B (i-P0)Er=i/ c 


d s 


0 [ s+e -(V5«+^ i )][i + ( i-p 0 )^( s ,O/C] 


(In 2)Au [(P on + SPt) +p 0 P off /(l -p 0 )] 


(34) 

The asymptotic characteristics of Q u are different from those 
of Qc . In heavy cell load, we have the scaling law of the 
“least” green user throughput, i.e., 0 (Ab/Au(-P O ii + <5Pt)) as 
p 0 -I- 0, which is dominated by the power consumption of the 
active mode and cell load. From the green user throughput 
point of view, heavy cell load should be avoided since any 
additional power consumption in the active mode completely 
gives rise to user throughput reduction. In light cell load, we 
have Qu = ©(Au/AbPos) as p 0 -+ 1, which is the same as 
the scaling law of the greenest cell throughput and also suffers 
the low throughput problem while achieving high greenness. 
Light cell load should also be dodged from the perspective of 
user’s quality of service. 

Since too heavy and too light cell loads do not benefit 
the green user throughput, we resort to finding the optimal 
cell load for Q\j if possible. Indeed, the following proposition 
shows that Q u also has a sole maximizer of cell load and how 
to find it. 


Proposition 6. For a given user intensity Au, there exists a 
/3 > 1 such that the optimal cell load that maximizes Q u in 
( |32| l can he accurately found by the sole fixed point of the 
following real-valued function 


ILu('u) 


=ttCAu 


(Pon - Ppff) 

<5P min r(i + 1 ) 



2 

ot + p 


(35) 


which means v* is the unique maximizer of cell load for Qy 
if and only (flLu^*) — v* = 0 and the optimal BS intensity 
is Ag = Au/u*. 


Proof: The proof is omitted here since it is similar to the 
proof of Proposition [5] ■ 

Obviously, the function in ( |35j ) for the optimal cell load of 
the green user throughput is similar to the function in ( |33j ) 
for the optimal cell load of the green cell throughput. Its 
fixed point also becomes larger if cell association schemes 
have a larger (. Under the same conditions, since green cell 
throughput is essentially smaller than green cell throughput 
and the increasing rates of Qc and Q u w.r.t. the cell load are 
pretty much the same in the light cell load, the optimal cell 
load for the green user throughput tends to appear earlier than 
that for the optimal cell throughput. This phenomenon can be 
observed from Proposition [5] as well as the simulation results 
in the following subsection. 


C. Simulation Results 

This subsection presents some simulation results to verify 
whether or not the previous results and discussions regarding 



Fig. 5. Simulation results of green cell throughput for nearest BS association 
and maximum received power association under Rayleigh fading and log¬ 
normal shadowing. The user intensity is Ay = 370 users/km 2 and n = 6 is 
used for eq. (31). 



— Analytical Sir in (33) (No Shadowing) 
Simulated Su (No Shadowing) 

— Analytical Su in (33) (Shadowing, <7, = 3dB) 
—B— Simulated Su (Shadowing, al = 3dB) 

■ Analytical Su in (33) (Shadowing, rr\ = 6dB) 
—^—Simulated Su (Shadowing, = 3dB) 

-Analytical Su in (33) (No Shadowing) 

- 4 - -Simulated Su (No Shadowing) 

-Analytical Su in (33) (Shadowing, o, = 3dB) 

-#■ -Simulated Su (Shadowing, <7, = 3dB) 

-Analytical Su in (33) (Shadowing, a* = 6dB) 

-A-Simulated Su (Shadowing, u, = 6dB) 


Xu/Xb (Cell Load) 


Fig. 6. Simulation results of green user throughput for nearest BS association 
and maximum received power association under Rayleigh fading and log¬ 
normal shadowing. The user intensity is Ay = 370 users/km 2 and n = 6 is 
used for eq. (33). 


the green cell and user throughputs are correct. We adopt 
the urban micro path loss model for 2GHz carrier frequency 
(140.7 + 37.6 log 10 d, d in km) given in |28) incorporated with 
average penetration loss (20 dB for 80% users that are indoor 
and 0 dB for the rest that are outdoor), antenna gain 5 dBi 
129]. The power consumption model of the small cell BS 
is based on the model of picocells in p4) . We summarize 
the common network parameters for all simulation figures as 
follows: P on = 6.8W, P off = 4.3W, P min = -106 dBm, 
S = 4.0, p 3 = 0 dB and a = 3.76. For simplicity, only 
two GCA schemes, i.e., nearest BS association and maximum 
received power association, are used to do simulation and 
compare the simulation results. 

Figs. 0 and [6] show the simulation results of the green 
cell and user throughputs, respectively. Both figures reveal 
that maximum received power association outperforms nearest 
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Fig. 7. Simulation results of the optimal BS intensity for green cell 
throughput under different user intensities. 



with the numerical results in Figs. [7] and [8] The correctness 
of Propositions [5] and [6] is validated. 

V. Conclusion 

In this paper, we investigate the optimal cell load, average 
and green throughputs from the cell and user perspectives in 
a small cell network with cell voidness. The impacts from the 
void cell phenomenon are two-fold. Firstly, it indicates that the 
previous works on throughput or rate are no longer accurate 
especially in a small cell network with dense BS deployment 
since no interferences are generated by the void cells. Sec¬ 
ondly, it inspires the idea of green power control that makes 
the void BSs switch to the dormant mode for saving power. An 
accurate expression of the void cell probability is successfully 
derived for GCA. The average cell and user throughputs that 
reflect the impact of void cells are theoretically characterized 
and they can be effectively boosted if appropriate cell load and 
GCA schemes with strong channel power (such as fading and 
shadowing) are adopted. The proposed green power control 
scheme not only exploits power saving on the void BSs 
but also gives birth to the idea of defining the green cell 
and user throughputs on the basis of the average cell and 
user throughputs. The green cell and user throughputs can 
be maximized by their unique optimal cell load that can be 
significantly improved by green power control and GCA. 

Appendix 

A. Proof of Theorem [7] 

The following proof essentially generalizes the proof of 
Lemma 1 in (23| . First, we define the void probability of 
<J>'nA for any bounded Borel set A' C M' 1 ’ as the probability 
that <&'(A) is equal to zero where T'fA) denotes the number 
of nodes in enclosed in A!. That is, 


Fig. 8. Simulation results of the optimal BS intensity for green user 
throughput under different user intensities. 


P [<b'(A) = 0] = exp 


A'(dX') 


BS association, shadowing power helps improve the green 
cell and user throughputs, and only one optimal cell load 
exists for every different simulation cases, as expected. Also, 
stronger shadowing power lets the optimal cell load happen 
later. To further explain this point, the simulation results of 
the optimal BS intensities for green cell and user throughputs 
versus different user intensities are provided in Figs. □ and [8] 
respectively. From these two figures, we see that the slopes of 
the green cell and user throughputs are initially both positive 
and reduce largely as the shadowing power increases (see 
the previous discussion for Qc in light cell load), which is 
because large shadowing power still has a strong influence on 
the coverage probability in ( |2T| i but has very less impact on 
the void cell probability. Propositions [5] and [6] also indicate this 
phenomenon and thus stronger shadowing power indeed leads 
to a smaller optimal BS intensity and larger green cell and user 
throughputs for the same user intensity, which is desirable. 
Finally, we observe that optimal BS intensities found by ( |33j ) 
with P £ (6,8) and ( |35| with /3 € (1,4) perfectly coincide 


because is a PPP. Let A, = Tj(A) where A C R d is a 
bounded Borel set and the probability that node A in $ is 
enclosed in Ai is 


Xi £ A i 


P [T i(X<) £ A] = P [TjTiiXD £ T?(A)] 

p[x'g(t fr i )- 1 (^)], 


where (TfT t ) _1 (A) = f A ,(TjT i )~ 1 {dX'). By definition, 
A(.4) is the average of <f?(4) and it can be calculated by 
applying the Campbell theorem as shown in the following: 


AM) 


[ A(dX) = E 

V i (Xi £ $ n a) 

Ia 

z —' \ / 

_Xi£$ 


= E 

= E 


E i^e^n^TO-'M 7 )) 


/(TTT)-i(A') 


P \X' £ $'] I'd' (dA) 
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= E 


= E 


1 


_ \J det(T' i T) 
1 

\J det(T T T) 


[ P[A" e &nA']v d (dX') 

J R d ' 

[ A'(dX') 


where (*) follows from the Jacobian determinant of two 
volumes. Thus, A(^4) = E . 1 = A '{A') and it follows 


that 

P[$'(-0 = 0] = exp I -E 


yj det(T T T) 
1 


\J det(T T T) 


)A 


A(dA’) 


for all Ci’s since nf=o C ^ represents the event that BS Bi 
is not associated by any users in its cell Ci and <l»u(C/) is the 
random number of users in cell Ci. The events £ij’s could be 
correlated each other such that the probability p@ is difficult 
to be explicitly computed. Instead, we find the bounds of 
because they can be characterized in closed-form as shown in 
the following. 

First, the lower bound on p@ can be found as follows. By 
considering a given and any two events £ij and £/,; are 

highly correlated under A t PI Aj ^ 0 such that P [£fj (T £j\ > 
P[£z.,]P[£;.j], we have the following 


= o 

P 

'-Su(Ci) 

n 

^u(Cz) 

= P 

£z,o n £;,i n • • ■ n £i ^, u ^ 

/). Since the void probability of a point 

_ j =0 





process completely characterizes the statistics of the process, 
$ is also a PPP. 

If d>' is homogeneous, then its intensity measure is 
A '{A!) = AV<j' ( A ') and thus the result in above can be further 
simplified to A(_4) = A'E y / det((T T T) _1 ) Vd(A) and it 

does not depend on the location of any node in $. As a result, 
$ is a homogeneous PPP with the intensity given in ([7]). 

B. Proof of Theorem [2] 

Suppose each user in cell Ci C R 2 of BS Bi has its own cell 
association region which is also a bounded Borel set on R 2 . 
For example, the cell association region of user Uj is denoted 
by Aj and all cell association regions contain BS Bi. Since 
all users in Ci adopt the GCA scheme, the probability that BS 
Bi is not associated by user Uj € Ci can be expressed as 


P [£ij] = lim 


sup {WtHiWUj- 


’ [£l,j\ = , ]im I 

v(Aj)—>00 


sup ha 

LB 4 €* B rUj\Bi 
> \Wj-Bi\ 


i - Bi 


= P Dj> D* 


-’Uj 


Dj is the distance from user Uj 


Pd> = E 


4>u(Ci) 

n 


j =0 


(Ci) 


$u(Ci) 

> n 

j=0 

Since the number of users in Ci is a Poisson random 
variable with parameter X\jv(Ci), i.e., P[$u(C;) = n] = 
{\ u v(Ci)) e -A U i/(C;) p or a gj ven the lower bound on p$ 

can be given by 

e ; n'Vir.ji )=e ( (3« 


3=0 


n— 0 \j =0 


gAu^(Cz) j-^\ 


According to the Slivnyak theorem 1221, the statistic property 


v(Aj)^oo |_B 4 e* B n.A J -\Bi 

> W t Hi\\Uj - Bi 


where £ 7 denotes the event that user Uj is not associated with 
BS Bi as its association region goes to infinity and v(Aj) 
needs to go to infinity in order to make P[£;j] include all 
BS association possibilities of user Uj. According to Theorem 
[I] and Remark 1, $b can be statistically mapped to another 
homogeneous PPP <1>b of intensity Ab = AbE[(JFFJ)S] since 
A' = Ab and T = (JF7T) _1 / Q in this case, and <1>b(-4j) stands 
for the number of BSs in (Aj (T <1>b) that is a Poisson random 
variable with parameter Xsi'iAj). Thus, it follows that 


evaluated at any point in a homogeneous PPP is the same such 
that P[£;j] is the same as P[£/,o] evaluated at user Uq (the 
origin) for all j. Accordingly, n; = oP[^-] in < |36l > is equal to 
(P[£i ;0 ])” and it follows that 

(Author 1 


n\ 


Pd, > E i (P[4ol)' 

l n =0 

= E |exp - P[£j, 0 ])j | • (37) 


in which Dq is the distance from 


and P [£i fi ] = P D 0 > D* 0 

the origin to the nearest BS in <I>b. The pdf of Dq is (x) = 
27t£Ab:e exp(—7 t£Ab 3; 2 ) 0301, whereby P[£;,o] is simplified as 


(38) 


P[£t,o] = 1 - E [exp (-<A B (A>) 2 ) 


Since all users in C / are uniformly distributed, we have 
fjj 0 ( r ) = and thus we can calculate P[£z,o] in closed- 
form as 


where Uj = ( WiHi ) 
to BS Bi, D* is the distance from Uj to its nearest BS in <I>b- 
Letting Ci be the cell of BS If . the void probability p@ can 
be found by calculating the following 


r'J'^U 2nrC 

P[£/,o] = 1 ^ / exp (-7rCA B r 2 ) :L ^- 1 dr 


v (Ci) 


= 1 - 


1 - exp(— (X B u(Ci)) 
Ab v(Ci) 


Clearly, 


> 1 — P[£;,o] and substituting this into ( |37j ) 


AB^fCi) _ 

results in the lower bound in (|T7J). 
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Next, the upper bound on p@ can be obtained as follows. 
The upper bound on P[£;j] is given by 


P[£;j] < lim 


v(Aj)—¥oo 


sup WiHiWUj-Bi 
B i £<S> B nAj\B l 

> WiHi inf {\\Uj — Bi\\~ c 

U j GG i 


*} 


'$u (Cl) 





n j 

$u (Ci) 

< (P 

Ri >D* 0 

) 

_ 3=0 






( b ) 
< 


(l - exp (-ttA b E (WiHi) ° f?o)) 
Ay 


E[4>u(Ci)] 


$u(Ci)/C 


1 - 


Ay+E[(WH)-.]A B J 

\ -<fu(Ci)/C 


= 1 + 


Ai 


1 + 


Au 


-CMj(Ci) 


D £ i 

3=0 


$u(G) 


< 1 + 


Ay 

CAb 


-C 


where Bq is the nearest BS in <J>b to the origin. Hence, for a 
given Wo Ho, ||Hq|| can be viewed as the distance from the 
nearest BS in to the origin that is scaled by 1/G. That 
means I 0 ||S 0 *II“ = (G-“/ 0 )(G“||H 0 *||“) ^ G-“I 0 ||H 0 ||“ if 
B 0 is the nearest BS in ‘I> B to the origin. For a given ||H 0 || = u 
and G, the Laplace functional of G~ a Iou a can be found by 


which is acquired by considering the farthest user of BS B[ in 
cell Ci. According to Slivnyak’s theorem, for a given <f>y(C;) 
we can have 


I 0 G-“(s|«) = E 
• sC " >; 


o -s(u/G) a I 0 


= E 


( 6 ) 

= e 


-E 


Jo” 1- 


B V. H z [u 2 /(\\Bi-i\\ 2 +u 2 )}2 
( uZ ) I 27rABrdr 


-3G~ a VHl 1+i 


where Ri denotes the distance from BS Bi to its farthest user 
in C,. Letting $y(C z ) 4 {Uj e C,. : (IF,//,) e C,} 

and Ho.j be the distance from BS Bi to the yth nearest user 
in Ci, we have 

’ =P[(H;) 2 > (D* 0 f 

iu(co 

E R h ^ ^0? 

3 =1 


WHiT^Ki ^ W 


(c) -7rABii 2 (l-P0)E|/“ ( 1- 


-sG~ a Hx 2 


drc 


where (a) is due to the fact that {||Hi|| 2 ,i = 1,2,...} are an 
one-dimensional Poisson process of intensity 7 tAb, (b) follows 
from the probability generation function of a homogeneous 
PPP 1221 and Theorem [T] and (c) is due to assuming the 
associated BSs are a thinning PPP of intensity Ab(1 — pf). 
By assuming Y is an exponential random variable with unit 
mean and variance, we have the following identity: 


f( 


1 - e~ sG aHx 2 )dx = 


sH 


r°° m 

sH 

/ p 

Y < 

Jo 

~ G a x 2 J 


d.T 


X - [ G a Y ) | dX Ey i(g q F 


V E[(WH)~*]\ b J 

( X y_\"* u(C,)/ «(c) 

V CAb7 - ' CA B/ 

where (a) follows from that fact that {( Ri) 2 ,l = 
1, 2, • • ■ } is an one-dimensional Poisson process of intensity 
AyE[(WH)°] and thus (Ri) 2 is equal to the sum of <f>y(C;) 
i.i.d. R 2 [31 j, ( b ) is due to first applying Jensen’s 

inequality on the term (WH) _ ° and then calculating the 
expectation regarding J?g, and (c) is obtained from (l + x) n < 
(1 + a;) 1 /" for x > 0 and n > 1 and E[$y(C;)] = 4>y(C;)/C > 
1. Therefore, we have 

'4>u(Ci) 


= (s7?) " G r 

which leads to the following result 

C u o, G -c, Io (s\u) = e — a b u 2 (i- P 0 )G- a /(«,G)_ (39) 

Then £|| B *||c«/ 0 (s) in ( |2Tj i can be obtained by calculating the 
following integral 


/*oo 

/ C uaG - aIo {s\u)e~ nXBu 2 d(TrX B u 2 ) 
Jo 


G 2 ^ 

G 2 + (1^7 0 )f(s,G 2 ) 


for any given ‘I'ljfC/) and the upper bound in 0 follows. 
This completes the proof. 

C. Proof of Proposition [2] 

According to (|8j, we can have the following CDF of ||Uq|| 

P[||H*|| < x] = P [||H 0 || < x(W 0 H 0 )-° 

= 1 — E [exp (—7 tAbG 2 ;e 2 ) , 


^||B 0 *|h/ 0 ( s ) = E 
= E 

Since £j|B*||a/ 0 (s) is concave for G 2 , we can use Jensen’s 
inequality to obtain sup{£|| B .||« /o (s)} = 1+(1 _ p0 W c)/c - 
Now we find 7c- According to ( p~3j ). Proposition IT] and the 
assumption that all associated BSs are a homogeneous PPP, 
the pdf of the area of associated cell Gg can be written as 

f«C5)(x) = (PAe( r ^ P f X)P e- pXB ^ )x (40) 
and the mean of \/v(Cq) can be found by 

E[1/KG 0 *)] = r -U(c 5 )(x)dx = pAe(1 ~/ 0) . (41) 

J 0+ x P — 1 













































14 


since p > Using this result, the average cell throughput in 
(|22l> can be expressed as 


Tc = E 


log 2 1 + 


Hn 


(ln2)(p- 1) 
Ab(i - ph) 


io\m 

Qt 

/»oo 


/ p 

Jo 

.11 

poo 


P) Jo 

P 


■E 


1 


//< 


o 


?5ll“Jo 
g 0 

ll^oll^O 


v(C* 

> (2 X - 1) 


da; 


> y 


Ay 


(V + 0' 


Moreover, we know 


P[£T > h] = 


—he z — 


(z-»sr 


V 2 n J~ 


dz 


= E 


— he 


= C e -z(h ), 


(42) 


Integrating both sides of the equation above and applying the 
appropriate boundary conditions of pi under the constraint that 
Pl is not a function of £(s, a) lead tcj^J 



=>w 2 = 


(Pon~ 
. (Pan 


Poff ) — tvV 2 -p (Ton 



-Poff) 


Since pi = 1 — (1 + v/p)~ p , ILc('y) in ( |33) > has a fixed point 
at v*. 

Now wee need to show the uniqueness of the solution of 
the differential equation in ( |45j ), i.e., the fixed point of Ii,c 
is unique. First, we rewrite ( |45) > as 

_dp®__ __ 

d [vf (1 +P0^(S,C)/C)] P oS - ( P on - Pos)p 2 /{s,C)IC 


where Z is a normal random variable with mean p s and 
variance of Thus, 7u can be expressed as 


Tc 


pA B (l-p 0 ) 

(ln2)(p-l) 


Ab(1 - pm) 

(In 2)(1 — 1/ p) 


r exp(-yI 0 \\B* 0 \\e- z ) 

L -FTP- dy 



C Io\\B*\\{s) 

(s + e~ z ) 


Using the above result and the definition of 7c- the average 
cell throughput can be further written as 


Tc = 


Ab( 1 -P 0 ) r Et^ + e- 2 )- 1 ] 


(In 2)(1 — 1 / p) 7 0 1 + (1-p 0 )^(s,C)/C 

Also, for any s > 0 we have 

pOO 


ds. (43) 


E 


1 

s + e 


-z 


1 


= £ 


/ °° / \ —1 

(s + e" (V e -P dx 

-OO ' ' 

(44) 


UJi 




where (d) is obtained by using a Gauss-Hermite quadrature 
integration with a sufficiently large integer n. Substituting ( |44| ) 
into ( [43] ) results in ( |22| ). For the case of no shadowing, letting 
Z in (|43[> be zero gives ([23]). This completes the proof. 


D. Proof of Proposition [5] 

The terms in Qc that pertain to the cell load can be lumped 
into the following sole function 


2J(u) = 


P0 


(1 +F 0 ^(s,C)/C)(F0(-Pon +KU 2 ) + (1 -_p 0 )P o ff) 


where v represents cell load, = 1 — p®{v) and k = 
(5P m i n F(l + -)/(7 tAuC)^ is a constant for a given Au- Then 
taking the derivative of Q3(u) w.r.t. v and making it equal to 
zero yield the following differential equation: 


^-7( S ,C)(Pon-Poff) 
- P® 


dpi 

du 


k£(s, C)v 2 d pi 


c 


dv 


1 + 


?( sX)P® 
C 


= [v* (1 +p 0 ('(s,C)/O] ■ 

(45) 


Assume both Ui (v) and ^(v) are two arbitrary solutions of 
© and let ip(v) = tpi — ip 2 - Since 0 < p 0 < 1, the following 
inequality 

n\ijj\dT 

I Poff — t(s, a)(P on - Poff) | 

always holds. By the Gronwall-Bellman lemma [32] , we know 
\tp\ < 0 and thus p = 0. Consequently, the solution is unique 
since any two arbitrary solutions tpi and are equal. 
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